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Abstract 


Searching for exoplanets with different methods has always been the focus of astronomers over the past few years. 
Among multiple planet detection techniques, astrometry stands out for its capability to accurately determine the 
orbital parameters of exoplanets. In this study, we examine the likelihood of extraterrestrial intelligent civilizations 
detecting planets in our solar system using the astrometry method. By conducting injection-recovery simulations, 
we investigate the detectability of the four giant planets in our solar system under different observing baselines and 
observational errors. Our findings indicate that extraterrestrial intelligence could detect and characterize all four 
giant planets, provided they are observed for a minimum of 90 yr with signal-noise ratios exceeding 1. For 
individual planets such as Jupiter, Saturn, and Neptune, a baseline that surpasses half of their orbital periods is 
necessary for detection. However, Uranus requires longer observing baselines since its orbital period is roughly 
half of that of Neptune. If the astrometry precision is equal to or better than 10 uas, all 8707 stars located within 
30 pc of our solar system possess the potential to detect the four giant planets within 100 yr. Additionally, our 
prediction suggests that over 300 stars positioned within 10 pc from our solar system could detect our Earth if they 
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achieve an astrometry precision of 0.3 jas. 
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1. Introduction 


More than 5400 exoplanets have been detected and 
confirmed to date (exoplanets.nasa.gov, 2023 July). Earth- 
sized habitable-zone planets turn out to orbit about one out of 
ten stars (Dressing & Charbonneau 2013; Petigura et al. 2013), 
and the search for life outside the solar system has experienced 
substantial impetus. Whether a planet is habitable or not 
depends on how far it is from the central star and its 
composition (Kasting et al. 1993; Gómez-Leal et al. 2018). 
Nowadays, more than 60 planets has been found to be habitable 
(Jones et al. 2006; Lovis et al. 2006; Anglada-Escudé et al. 
2012; Barclay et al. 2013; Borucki et al. 2013; Tuomi et al. 
2013; Robertson et al. 2014; Jenkins et al. 2015; Dittmann et al. 
2017; Gilbert et al. 2023), most of which are detected by the 
transit and radial velocity method. Neither the transit nor radial 
velocity method provides complete physical parameters of one 
planet, and both methods prefer to detect planets close to the 
central star. On the contrary, the astrometry method can 
provide three dimentional characterization of the orbit of one 
planet (Perryman et al. 2014; Wu et al. 2016) and has the 
advantage to detect planets far away from the host star. 

To date, only one giant planet has been detected by the 
astrometry method (Sahlmann et al. 2013) because of the 
limitation of detection precision. The detection of a habitable 
Earth-sized planet orbiting around a Sun-like star located 10 pc 
away from us would require a precision of sub-jias, which is 
hardly achieved by the current astrometry observation such as 


Gaia (Perryman et al. 2014). However, it is very promising in 
the near future with the coming of a new era with high 
astrometry precision of uas (Yu et al. 2019; Ji et al. 2022; Jin 
et al. 2022; Tan et al. 2022). 

Here we propose a probing question that supposing the 
extraterrestrial observers are using the astrometry method and 
are also surveying the Galaxy for habitable worlds, which of 
them could discover the planets in the solar system, even the 
Earth? Previous works have investigated the region in which 
the Earth will be observed transiting in front of the Sun (Heller 
& Pudritz 2016; Kaltenegger & Pepper 2020; Kaltenegger & 
Faherty 2021) and the frequency that the Earth will be detected 
by other civilisations through photometric microlensing 
(Suphapolthaworn et al. 2022). 

In this work, we study the possibility that extraterrestrial 
lives detect the planets in the solar system via the astrometry 
method with different observational precisions. We describe 
how we simulate astrometric data in Section 2. In Section 3, we 
present how to identify planetary signals and how to fit the 
orbital parameters of the planets. The detection of the four 
giants in the solar system by nearby stars are discussed in 
Section 4. We briefly conclude our results in Section 5. 


2. Simulation of Astrometric Data 


The astrometry method measures the movements of the stars 
projected onto the celestial sky. Following the method 
described in previous works (Black &  Scargle 1982; 
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Wu et al. 2016; Yu et al. 2019), the projected movement of the 
star in the R.A. (x) and decl. (y) at time t can be modeled as: 


x(t) = xo + H(t — to) — Pow + X(t) + ox (1) 


and 


y(t) = yo + ny — to) — Er + Y(t) + oy, (2) 


where xo and yo are the coordinate offsets, ux and py are the 
proper motions of the star, P, and P, are the parallax 
parameters which will be provided by the observation. 7 is 
the annual parallax of the star. X(f) and Y(t) are the movements 
of the host star around the barycenter of the system due to the 
planetary companions. c, and o, are single-measurement 
astrometric errors. 

In our fiducial simulations, we made several assumptions 
regarding the extraterrestrial observer's location and observa- 
tional parameters. As part of our simulations, we placed the 
observer at a distance of 10 pc from our Sun. The observer 
orbits its central star in a circular orbit with a period of 1.25 yr 
and measures the coordinates of the Sun (x(t) and y(t)) every 
0.2 yr. We also conducted simulations with data cadence of 
0.1 yr and found that the results change very little. To account 
for the astrometry precision, we assumed that the observer has 
a measurement uncertainty of 10 jas. Therefore, the individual 
coordinate uncertainties c, and o,, were chosen from a 
Gaussian distribution with a median value of 0 and a standard 
deviation of 10 yas. For the coordinate offsets xo and yo, we 
assume both to be 10 mas. Additionally, the proper motion of 
the Sun with respect to the observer is assumed to be 
50masyr ! and -30masyr ! for the x and y directions, 
respectively. To model the parallax effect, we used the 
observer’s orbit and defined the functions P, and P, as follows: 
P(t) = sin(l.6zt + 0), P(t) = cos(l.6zt + ¢), where @ 
represents the orbital phase of the observer. Finally, we 
assumed an observing baseline of 170 yr for the simulations. 

The movements of the Sun due to the presence of the eight 
planets X(t) and Y(t) are simulated using the REBOUND code 
(Rein & Liu 2012). All eight planets in the solar system are 
included. The orbital parameters of the planets are given by the 
JPL Solar System Dynamics website, with respect to the mean 
ecliptic and equinox of J2000. In our fiducial simulation, the 
line of sight of the extraterrestrial observer is assumed to be 
perpendicular to the mean ecliptic. We integrate the solar 
system over a duration of 170 yr and record the coordinates of 
the Sun (X(t) and Y(f)) relative to the barycenter of the solar 
system every 0.2 yr. 


3. Planetary Signal Identification and Orbital Para- 
meter Fitting 


Assuming that we are extraterrestrial civilizations and we 
have measured the movement of the Sun for 170 yr. Now we 
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analyze the data to see if we have any detection. Although we 
have included the gravitational interaction between planets 
when we simulate the astrometric data of the host star, it is 
ignored when we fit the orbital parameters since it has little 
influence on the motion of the host star (Sozzetti et al. 2001; 
Casertano et al. 2008). In our parameter fitting procedure, X(t) 
and Y(t) are modeled as (Catanzarite 2010): 
i=N 
X(t) = > (cos E; — e)A; + 


i=1 


1 — e? (sinE)F, (3) 
and 
i=N 
Y(t) = Y (cos E; — e)B; + 41 — e? (sinE)Gi, (4) 
i=j 


where N is the number of planets orbiting around the central 
star, i represents the in planet, E; is the eccentric anomaly, e; is 
the orbital eccentricity, A; F;, B; and G; are Thiele-Innes 
constants, given as: 


A; = a;(cos Q; cos w; — sin Q; sin w; cos [;), 

Bj = a;(sinQ; cos w; + cos Q; sin w; cos £j), 

F; = a;(—cos Q; sin w; — sin Q; cos w; cos Å), 

G; = a;(—sin Q; sin w; + cos Q; cos w; cos I), (5) 


where a; is the astrometric signature of the host star due to the 
reflex motion in the presence of the in planet. Q, w and J are the 
longitude of ascending node, arguments of pericenter and the 
orbital inclination of the planets, respectively. 

We search for the planetary signal and then fit the orbital 
parameters of the planets following the procedures as we 
described in Wu et al. (2016). Here we briefly describe the 
steps. 

Step 1, ignore the planetary influence on the star and use the 
linear least squares method to fit the five stellar parameters xo, 
Yo. Hx» Hy and T. 

Step 2, remove the coordinate offsets, stellar proper motion 
and parallax from the data and search for periodical signals in 
the residuals using the Lomb-Scargle periodogram (Black & 
Scargle 1982). We calculate the periodogram of the residuals in 
the x and y directions and record the most significant peak in 
each direction. Then we choose the peak with smaller false 
alarm probability (FAP). If the peak has an FAP < 10 ^, we 
claim that we have identified a certain planet signal and the 
corresponding orbital period is adopted as P4. 

Step 3, fit xo, Yo, Hx» My, T, P1, €1 and fo}. toj is the perihelion 
moment of the planet. This is processed via the Levenberg- 
Marquardt (LM) algorithm (Marquardt 1963) and the Markov 
Chain Monte Carlo (MCMC) fitting procedure. Given Pi, e; 


and f9;,, the terms cos =Æ — e; and 4/1 — e (sin Ej) can be 
determined. Then Equations (3) and (4) are easily inverted by 
linear least squares to yield the Thiele-Innes constants. The 
motions of the Sun produced by the planet are calculated using 
Equations (3) and (4). Together with the five stellar parameters, 
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we can calculate the fitted projected motion of the star using 
Equations (1) and (2). We first fit xo, Yo, fx, Hy, T, P1, e, and toi 
using the LM method. Initial values of xo, Yo, Hx» Hy, and Pı 
are given by Step 1 and 2, while e; is randomly chosen between 
O and 1, t9; is randomly chosen between 0 and P,. The LM 
fitting process is repeated for 100 times. Then we choose the 
best-fit parameters with the smallest reduced x? as initial values 
of the following MCMC fitting procedure. We adopt the open- 
source Python package emcee (Goodman & Weare 2010; 
Foreman-Mackey et al. 2013) to sample the parameter space 
and estimate the posterior distribution of parameters. We run 
emcee with 64 walkers for 20000 + 30000 x N iterations (N is 
the number of planets identified). The initial positions of the 
walkers are drawn from Gaussian distributions with median 
values given as the best-fit parameters of the LM fitting process 
and standard deviations of 10 ? to accelerate the fitting 
process. We conduct autocorrelation analysis and find that all 
chains are converged in our fitting procedure. 

Step 4, remove the coordinate offsets, proper motion, 
parallax and stellar motion due to the planet identified in Step 
2 using the best-fit parameters calculated in Step 3 from the 
original astrometric data. Then we continue to search for 
periodic signals in the new residuals. If there is one, then we fit 
the data with a two-planet reflex motion model. 

Step 5, repeat Step 2 to Step 4 until no periodic signals are 
identified. 

In our fitting procedure, we have a total of 5+3xWN 
parameters to fit since we have assumed Keplerian orbits for 
each planet, which largely reduce the parameters to be fitted 
and ensure the parameter precision at the same time. The 
semimajor axis of the planets can be obtained using the 
Kepler’s third law giving the orbital period of the planets, while 
the planetary masses are calculated via m;a; = modo, where 
m, is the solar mass (assumed to be precisely determined via 
other methods by the extraterrestrial intelligence, like the 
spectrometry or astroseismology) and a is the semimajor axis 
of the Sun when orbiting around the barycenter determined by 
the Sun and the ij, planet (which is obtained in Step 3). Readers 
are refereed to Catanzarite (2010) and Wu et al. (2016) for 
more detail. 

In Figure 1, we present the data residuals and power 
spectrum for each step of the fitting process in our fiducial 
simulations. After Step 1 and Step 2, the power spectrum of the 
data residuals (labeled as 1,, O — C) exhibits a prominent peak 
at approximately 11.9 yr, indicating the successful identifica- 
tion of Jupiter. Then we move to Step 3 and Step 4, the updated 
data residuals (labeled as 2.4 O — C) are also shown in Figure 
1, with their power spectrum peaking around 28.8 yr. This peak 
signifies the detection of Saturn. Continuing this iterative 
process, we repeated the aforementioned steps, ultimately 
leading to the identification of Neptune and Uranus. After the 
detection of the four giants, no peak with FAP < 10^ appears 


in the final residuals, suggesting that none of the small planets 
in the solar system is detectable in our fiducial simulations. 


4. Results 


4.1. The Characterization of the Four Giants in the Solar 
System 


The amplitude of the astrometric motion of the Sun produced 
by a planet with mass m, and semimajor axis a observed by an 
observer with a distance of d is: 


EP ELO EA CR 
EE E29 Fed Eod and "d 


With an observing baseline of 170 yr and an observational 
error down to 10 as, all the four giants are successfully 
detected and characterized. This is expectable since the signal- 
noise-ratios (S/Ns, defined as a/o, where o is the 
observational error) of the four giants calculated using 
Equation (6) are far larger than 3, according to the detection 
criterion given by Wu et al. (2016). Other small planets in the 
solar system are hardly detectable because of their small S/Ns. 
We show the posterior distributions for all parameters that are 
fitted in Figure 2. The first half iterations are thrown as burn-in. 
We find that the five stellar parameters (xo, yo. Hx» Hy and 7) and 
orbital parameters of the inner three giants (Jupiter, Saturn and 
Uranus) are all well-constrained with nearly Gaussian 
distributed posteriors, indicating that the parameters converge 
well. For the outermost planet Neptune, the orbital eccentricity 
(e3) and perihelion moment (fo3) are not well-constrained since 
the planet only finishes one complete orbit during 170 yr. The 
planetary mass and semimajor axis of the planets can be easily 
calculated using the best-fit parameters as we have mentioned 
in Section 3. All the four giants are well characterized with 
relative fitting errors smaller than 1% for both the orbital period 
and planet mass. The relative fitting error of parameter 0 is 
given as €g = [0g — Orucl/Ptruc- Ont is calculated as the median 
value of the posterior distribution of 0. 


4.2. The Influence of Observing Baseline and Observa- 
tional Error 


We also investigated the detection of the four giants with 
different observational errors and observing baselines. In our 
fiducial simulations described in Section 2, we fix the 
observing baseline to be 170 yr and the observational error to 
be 10 uas. Now we gradually decrease the observing baseline 
from 170 to 10 yr, with a step of 20 yr. To account for the 
detection of our solar system by missions like Gaia (Perryman 
et al. 2014) and CHES (Ji et al. 2022), we further extend the 
observing baseline down to 4 yr. We also considered different 
observational errors: 1, 3, 10, 30, 100, 300, 1000, 3000 and 
10,000 uas. Other assumptions such as the distance of the 
observer and the sampling cadence remain the same. In the new 
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Figure 1. Data residuals and power spectrum of the four giants after each fitting step. Top: the simulated astrometric data on the x (shown in red) and y (shown in blue) 
directions. Left: The data residuals after each fitting step. Right: the power spectrum of the data residuals shown on the left. 


simulations, we assume that the coordinates offsets, proper 
motion and parallax of our Sun are already well determined and 
carefully removed from the astrometric data by extraterrestrial 
intelligence before the fitting process starts. This assumption 
will largely reduce the computational time. (We have 
conducted a small group of simulations including the 
coordinate offsets, proper motion and parallax, and we find 
that they have little influence on the characterization of the 
planets.) Then we start the fitting procedure as described in 
Section 3 but now we skip step 1. 

We show the relative fitting errors of planet mass as a 
function of the observing baseline and observational error for 
each of the four giants in Figure 3. Only planets with relative 
fitting errors of orbital period smaller than 0.1 (ep < 0.1) are 
shown. For planets with ep<0.1, their planetary mass are 
mostly well-fitted with e,, < 0.1. Only a few exceptions with 
large observational errors or short observing baselines have 
large fitting errors of planet mass. There are several cases that 
planets are detected with ep 0.1, however, their planet mass 
are generally poorly fitted with e,, > 1. Therefore, we claim a 
planet is well characterized if it has ep « 0.1. 

As we have pointed out in Wu et al. (2016), the detection of 
a planet using astrometry method relys on the S/N of the planet 
and the observing baseline. Here we show the contours of the 


S/Ns in Figure 3. We find that all the four giants can be 
successfully detected and well-characterized as long as their 
S/Ns » 1 and the observing baseline exceeds 90 yr. In general, 
planets with S/Ns > 1 and observing baseline longer than 
about half an orbital period could be detected. However, this is 
not the case for Uranus. Because the fitting of the orbital period 
of Uranus is largely influenced by that of Neptune, whose 
orbital period is about two times that of Uranus. Not until the 
orbital period of Neptune is successfully identified will Uranus 
be well-characterized. There are exceptions that planets are 
detected with S/N « 1. However, these detections generally 
have large fitting errors of planet mass. 

If astrometric missions conducted by extraterrestrial civiliza- 
tions are similar to Gaia or CHES, with typical observing 
baselines ranging from 5 to 10 yr, the detectability of the giant 
planets in our simulations would be limited. Specifically, only 
Jupiter would be detectable under these circumstances. 
However, if Gaia were able to achieve an observational error 
down to 100 as, it would be possible to characterize Jupiter 
with an accuracy of e,, < 0.5. Alternatively, if CHES, with an 
observational error down to 1 pas and and an observational 
time of approximately 6 yr, conducted the mission, Jupiter 
could be characterized with an accuracy of €m ~ 0.3. 
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Figure 2. Posterior distribution of xo, yo, Hx Hy, T, P1, €1, to,1, P2, e2, to,2, P3, €3, to,3, P4, e4, and to4. Four planets are detected by extraterrestrial intelligence located 
10 pc away with an observing baseline of 170 yr and observational precision of 10 pas. 


4.3. Which Stars Could Detect the Four Giants in the identify 8707 stars from the Gaia Catalog of Nearby Stars (Gaia 
Solar System Collaboration et al. 2021) that lie within 30 pc to the solar 

We move forward to estimate how many neighboring stars in system. We calculate the S/N of each giant planet observed by 
the Galaxy could detect the four giants in our solar system. We each star with different assumed observational errors. We find 
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Figure 3. The relative fitting errors of the planet mass €, as a function of observing baseline and observational error. Different colors represent different c,,. Gray 
circles represent none detection of planet signals or the identified planet has large fitting errors of orbital period (i.e., e» > 0.1). The gray dashed lines represent the 


contours of different S/Ns. 


that all 8707 stars have the possibility to detect and well- 
characterize the four giants if they could achieve an astrometric 
error down to 10 pas and observe the solar system for enough 
long time (such as 90 yr). If the observational error is as large 
as 100 pas, only 183 neighboring stars could detect all the four 
giants, but all of them could detect the Jupiter within 10 yr. We 
also estimate the number of neighboring stars that could detect 
our Earth. About 310 neighboring stars located within 10 pc 
from our Sun have the potential to detect the Earth if the 
observational error is as small as 0.3 uas. With a larger 
observational error, such as | pas, only 8 stars located within 3 
pc from the Sun could possibly detect the Earth. 


5. Conclusion 


In this paper, we study the possibility that extraterrestrial 
intelligence detect the planets in our solar system. We find that 
all the four giants in our solar system could be detected and 
well-characterized as long as they are observed for at least 90 yr 
with S/N > 1. For all 8707 stars lying within 30 pc to the solar 
system, all of them have the potential to detect the four giants 
within 100 yr if they could achieve an observational precision 
down to 10 yas. If the astrometry method can achieve sub pas 


precision like 0.3 pas, then even our Earth will be detectable by 
extraterrestrial intelligence. 

In each of our simulations, we assume a constant 
observational error during the long observing baseline for 
simplicity. A more reasonable assumption may be a decreasing 
observational error with the increase of observing baseline. 
Besides, the sampling cadence is fixed to be 0.2 yr in our 
simulations. This is hardly achieved in real observations. 
Finally, our simulations is truncated at 170 yr since longer 
observing baseline requires longer computational time to fit the 
orbital parameters of the planets. However, we expect that 
longer observing baseline would allow the detection of plants 
with larger observational errors. These should be further 
considered in future works. 

This study primarily addresses the likelihood of extra- 
terrestrial civilizations in the vicinity of our solar system 
detecting our own system. However, it is important to note that 
the existence of extraterrestrial lives remains uncertain, and if 
they do exist in planetary systems similar to ours, their 
presence could be incredibly rare. Cumming et al. (2008) 
demonstrated that the occurrence rate of cold Jupiters around 
stars similar to the Sun is only 10%. Consequently, the chances 
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of our solar system being in proximity to a significant 
population of extraterrestrial civilizations are currently 
considered to be very small. 
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